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Abstract 

We study the implications of non-standard Higgs Yukawa couplings to light quarks on Higgs- 
portal dark matter phenomenology. Saturating the present experimental bounds on up-quark, 
down-quark, or strange-quark Yukawa couplings, the predicted direct dark matter detection scat¬ 
tering rate can increase by up to four orders of magnitude. The effect on the dark matter annihila¬ 
tion cross section, on the other hand, is subleading unless the dark matter is very light - a scenario 
that is already excluded by measurements of the Higgs invisible decay width. We investigate the 
expected size of corrections in nnilti-Higgs-doublet models with natural flavor conservation, the 
type-II two-Higgs-doublet model, the Giudice-Lebedev model of light quark masses, minimal flavor 
violation new physics models, Randall-Sundrum, and composite Higgs models. We find that an 
enhancement in the dark matter scattering rate of an order of magnitude is possible. Finally, we 
point out that a discovery of Higgs-portal dark matter could lead to interesting bounds on the 
light-quark Yukawa couplings. 
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I. INTRODUCTION 


In Higgs-portal models [1-10] of dark matter (DM) the Higgs is usually assumed to be 
completely Standard Model (SM) like apart from its couplings to DM. Experimentally, only 
the couplings of the Higgs to the heaviest particles of the SM are currently well constrained. 
The couplings to gauge bosons are found to be in agreement with the SM predictions at the 
(9(20%) level, while the constraints on the couplings to third-generation fermions are some¬ 
what weaker [11, 12], Much less is known experimentally about the couplings of the Higgs to 
the first two generations of fermions. The couplings to u, d, s, and c quarks could be as large 
as the SM bottom Yukawa coupling or be absent altogether [13-16]. The Higgs couplings 
to top and bottom quarks will be quite well known by the end of the high-luminosity LHC 
run. Some progress is also expected on the measurements of Higgs couplings to charm and 
strange quarks [13-15]. 

Large u-, d-, and s-quark Yukawa couplings, comparable in size to the 6-quark Yukawa, 
generically require fine-tuning. A large Yukawa coupling implies a large contribution to the 
quark mass from the Higgs vacuum expectation value (vev). This would then need to be 
cancelled by a different contribution to the «-, d-, and s-quark masses, unrelated to the 
Higgs vev. The opposite limit, where the observed Higgs does not couple to the light quarks 
at all is easier to entertain. It simply requires a separate source of the light-quark masses 
(for an extreme example see, e.g., [17]). 

Modified light-quark Yukawa couplings could, in principle, have important implications 
for DM phenomenology. In this article we investigate how the Higgs-portal DM predictions 
change if the Higgs couplings to the light quarks differ from the SM expectations. We 
first allow for an arbitrary flavor structure of the Higgs Yukawa couplings, only requiring 
that they satisfy the current experimental bounds. In Section II we derive the implications 
for direct DM detection, indirect DM detection and the collider searches. We show that 
vanishing couplings of the Higgs to light quarks only have a relatively small impact on 
these observables. Saturating the loose current bounds on the light-quark Yukawa couplings 
would, on the other hand, lead to drastically enhanced scattering cross sections on nuclei 
while leaving the relic density and annihilation cross sections nearly unmodified. 

Clearly, an enhancement of the light-quark Yukawas by factors of (9(100) or more, as 
allowed by current data, requires considerable fine tuning of the quark-mass terms and 
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hence seems quite unlikely. In Section III we, therefore, explore the deviations in the Higgs 
Yukawa couplings for a number of beyond-the-SM scenarios and flavor models. This leads to 
more realistic expectations as to how large the deviations in the direct DM detection rates 
can be due to the poorly known Higgs couplings to the light quarks. Note that we assume 
the DM to be a flavor singlet and that the new flavor structure of the interactions with the 
visible sector is only due to the modification of the SM Higgs couplings. DM that is in a 
nontrivial flavor multiplet has been investigated in [18-29], while our study is closer in spirit 
to the work in [30-32] where the flavor dependence of the DM signals for flavor-singlet DM 
has been explored. 

A somewhat surprising result of our investigation is that, if DM is discovered and turns 
out to be a thermal relic predominantly interacting through a Higgs portal, it could be used 
to constrain the light-quark Yukawa couplings. This is discussed in more detail in Section IV. 

We summarize our results in Section V. 


II. HIGGS PORTAL WITH NON-TRIVIAL FLAVOR STRUCTURE 


We assume that DM and the SM fields are the only light degrees of freedom. The 
remaining new physics (NP) particles can be integrated out so that one can use an Effective 
Field Theory (EFT) approach. The couplings of DM to the SM are given by the Higgs-portal 
Lagrangian 


9xX ] xH ] H , 


scalar DM; 

fermion DM; (1) 

vector DM. 

Above, the fermion DM can be either a Dirac or Majorana fermion (in either case we use 
four component notation). After electroweak symmetry breaking (EWSB) we have 


1 1 

C x= 9xjXXH ] H + ig x -xisXH ] H, 


H ] H 


1 

2 


(v'w + 2 v w h + h 2 ), 


( 2 ) 


where Vw = 246 GeV is the vacuum expectation value (vev) of the Higgs field. The above 
interactions therefore lead to annihilation of DM into both single Higgs, XX h, and double 
Higgs, XX hh, final states. 

The Higgs-portal operator for fermionic DM has mass dimension five and is suppressed 
by the new physics scale A. The Higgs-portal interaction for fermionic DM can also be 
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re-written as C x = (g x + ig x )xLXR^ H / A + h.c.. For g x ^ 0 the interaction is thus both 
P- and CP-violating. The interaction for vector DM is most probably also due to a higher¬ 
dimensional operator in the full theory. For instance, if Xu arises from a spontaneously 
broken gauge symmetry in the dark sector, then g x ~ v 2 D /A 2 , where vd is the vev of the field 
that breaks the dark sector gauge invariance, while A is the mass of the mediator between 
DM and the Higgs. 

The relevant terms, after EWSB, in the effective Lagrangian for the Higgs couplings to 
the SM particles are given by 



m e n, 


h 


v w ^ v w 


(3) 


where the k * are real. A sum over the SM quarks, q = u,d, s, c, b, t, and charged leptons, 
£ = e,n,T is implied, and we have assumed custodial symmetry. The h —* 77 coupling is 
not relevant for DM phenomenology, since its effects are suppressed compared to the Higgs 
couplings to gluons . 1 The couplings are normalized such that K q = Kg. = Ky = K\ = 1 
correspond to the SM. The experimental constraints on the couplings of the light quarks 
to the Higgs, obtained from a global fit to current data, are < 0.98 mb/m u , |/Cd| < 
0 .93mb/md, |/c a | < 0.70 m^jm^ where only one of the light Yukawa couplings was left to 
float in the fit, while all the other Higgs couplings are set to the SM values [13]. Higgs 
couplings to the light quarks of a size comparable to the coupling to the b quark are thus 
still allowed. In (3) we do not allow for flavor violating Higgs couplings, since these are 
already tightly constrained from both Higgs decays and low-energy observables [33-36]. 

In the SM, keeping the Higgs on shell, the hGG coupling arises predominantly from the 
one-loop top-quark contribution. The in (3) encodes only the potential NP contribu¬ 
tions, and vanishes in the SM. I 11 the global fits a parameter K g is introduced that gives the 
total h —» gg amplitude, including the SM contributions [11, 12]. We have (see, e.g., [37]) 


1.03/ct + < P 


(4) 


It could be relevant for direct detection if the scattering on electrons dominates. This requires very light 
DM, of order the electron mass. Such light Higgs-portal DM is excluded by the constraints on the Higgs 
invisible branching ratio. 
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Figure 1: DM annihilation channels in the Higgs-portal models. 

At present, significant CP-violating Higgs couplings to fermions and gluons are still 
allowed experimentally (see, e.g., [38]), so we also discuss their effect on the Higgs interactions 
with DM: 

Afgcpv = -ik^q^qh - ike—hdh + R* p -^hG%G a > JI '. (5) 

Here, the ki are real parameters; in the SM, we have ki = 0. Moreover, G a ^ u is the gluon 
field-strength tensor and G a ^ u = G° t6 its dual. The normalization of the hGG term is 

chosen such that integrating out the top at one loop one obtains k g = k t + k^ p . Accordingly, 
we have Br (h —)■ gg ) cx n 2 g + {2>k g /2) 2 . The CP-violating couplings of the Higgs to ZZ and 
WW are already well constrained, and we thus set them to zero. In the numerical analysis 
below, we will also assume k t — 0 , for simplicity. 

The modified Higgs couplings change the usual Higgs-portal predictions for DM annihi¬ 
lation rates, the relic abundance, and direct detection rates. In the following, we discuss 
these modifications in detail. 

A. Annihilation cross sections 

The dominant DM annihilation cross sections in the Higgs-portal models are XX ~^ 
bb,W + W ~, ZZ,tt, and XX hh. The first four proceed through the s-channel Higgs 
exchange, while XX hh receives additional contributions from t- and tt-channel y exchange 
as well as from the four-point contact interaction (cf. Fig. 1). The XX bb channel is 
only relevant if the other channels are not kinematically allowed, i.e. for light DM masses, 
m x < rn w . 

Tl ie XX bb annihilation cross section assuming SM Higgs couplings is given for scalar 
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(S), Dirac fermion ( DF ), and vector (V) DM by 


{ a bb V ^)sM 

( a ti V rel) SM 


N r 


almlPl 


47 T ( s ~m 2 h ) +m 2 h T 2 h 
Nc mI g 2 x (s - 4 m 2 x ) + ~g\s 
8n A 2 
,2 


N c 9 


s-ml) +m 2 h T 2 h 
2 ' 3 


A 3 


b ? 


X m b n3 

Pb 


9 167 T m% 

A. 


r 


x + 4 r x 


(s - ml) 2 + m 2 h r 2 ’ 


( 6 ) 

(7) 

( 8 ) 


where, here and below, yfs is the center-of-mass energy, = 4m|/ s, (3k = x/1 — rk is the 
velocity of particle k, and n re i = 2/3 x is the relative velocity of the DM particles. If the Higgs 
coupling to the 6 -quarks differs from the SM value, the annihilation cross section is rescaled 
as 


<Tbb= (*b + ii b/Pb) 0 'l 5 M - (9) 

The annihilation cross sections cjfj to the other fermions are obtained with the obvious 
replacement 6 —> f in the above expressions. Since the Higgs couplings to the light quarks 
are poorly constrained experimentally, the DM annihilation to two light quarks can be 
comparable to XX bb and can be important for light DM, m x < mw- 

For heavy DM, rn x > mw, the annihilation cross-sections into a pair of W or Z bosons 


are 


„ _ k .2 SM 

°VV — KyCTyv , 


( 10 ) 


V — W, Z. The annihilation cross sections assuming the SM Higgs couplings to W are given 
by 


(°WW' y rel)sM — V ( 1 ~ T W + ~T 


8tt 


w 


&WW V *d)sM — 


167tA 2 
,2 


(3 W ( 1 — r w + -r 2 


s - m 2 ) 2 + m 2 T 2 ’ 

s\g 2 x (s-4m 2 x )+g 2 x s\ 


w 


cr ww v rel)sM — 


9 X 


288n 


Pw ( 1 — fw + T r w 


(s - m 2 ) + m 2 T 2 


1 — r Y + -r( 


(s - ml f + m\Tl ' 


( 11 ) 

( 12 ) 

(13) 


for scalar, Dirac fermion, and vector DM, respectively. The \X —> ZZ annihilation cross 
section is obtained by replacing W —> Z, and multiplying all expressions by an extra factor 
of 1/2 since one has two indistinguishable particles in the final state. 

The XX hh annihilation cross sections for scalar, Dirac fermion, and vector DM, are 
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given by 


s Phg 

VhhVrel = 


2 r 

x 


64n ml 

X. *- 


1 + 


3 KoMl 
4m 2 - Ml 


2v w9x 

Ml - 2 m 2 


1 2 


DF V Ph{9 2 x + 9 2 xPx) 

a hh Del = 


(J hh V rel 


32vrA 2 
Ph 


1 + 


3aca Ml 


+ 


4 9x m x v w 


5767rm 2 


3<? 


4m 2 -Ml A(2m 2 -M 2 )_ 
3 n\Ml 


(14) 

(15) 


4m 2 - Ml 


+ 1 + 


3„,2 


lQg x v 


x u w 


2m 2 


3k\M, 


Ml \ 4 ml 


Ml 


4 g\ v w 

(2 ml - Ml) 

+ 1 


6 - 


4 Ml | Ml 

m x m x 


Ml\ 


4m 2 / 


(16) 


In this result we display only the leading terms in the expansion in powers of fi x . The 
contribution of the s-channel Higgs exchange diagram is proportional to the rescaling of the 
trilinear Higgs coupling, n\. The latter is completely unknown experimentally, at present, 
but can be measured to 0(20%) at the end of the LHC [39, 40]. Since we are mostly 
interested in the effects of flavor modifications we will set it to the SM value, K\ — 1, in the 
numerics below. All cross sections for Majorana DM can be obtained by multiplying the 
corresponding Dirac DM cross sections by a factor of 4. 


B. Relic abundance 

The DM relic abundance Ddm is proportional to l/<Tn re i, where a is the annihilation cross 
section. Assuming that the DM in our scenario accounts for all of the observed relic density, 
the measured value Q, DM h 2 = 0.1198(26) [41] fixes g x for a given value of m x . The resulting 
constraint in the m x - g x plane is denoted for the different cases by a red line in Figs. 2 to 9. 
In Fig. 2, we compare two limits of the Higgs portal for the scalar DM: the case where the 
Higgs does not couple to the light quarks at all (left panel) to the case where the Higgs has 
SM Yukawa couplings (right panel). The two relic abundance curves coincide apart from 
very light DM, with m x below the charm and tau threshold. If such light DM did not couple 
to the u, d and s quarks, this would result in noticeably reduced annihilation cross sections 
and, thus, in larger relic abundance. In both cases, the dominant annihilation process is 
still given by XX h* 99- For very light DM the correct relic abundance is obtained 
only if the coupling of the Higgs to DM, g x , is nonperturbatively large. The yellow regions 
in Fig. 2 denote the value of g x for which the total Higgs decay width would be larger than 
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Scalar DM, zero light Yukawas 


10 10 2 10 3 10 4 
m x [GeV] 


Scalar DM, SM Yukawas 


10 10 2 10 3 
m x [GeV] 


Figure 2: Bounds from LUX (blue band), XENONIOO (green band) and the invisible Higgs decay 
width (black dashed line and grey region) on the Higgs-portal coupling g x for scalar DM, assuming 
vanishing (left) and SM (right) Yukawa couplings to u, d, s quarks. The red line denotes g x as a 
function of DM mass, rn x . for which the correct relic abundance is obtained, while g x in the yellow 
region leads to non-perturbatively large Higgs decay width, T/, > rrih, and is excluded. Constraints 
from Fermi-LAT searches for DM annihilation in dwarf spheroidal galaxies are denoted by the 
orange band. 

its mass, T h > rrih, and are thus excluded. 

The same comments apply to the case of vector DM, shown in Fig. 4, CP-conserving 
Dirac fermion DM, shown in Fig. 6, CP-violating Dirac fermion DM, shown in Fig. 8, and 
also for Majorana fermion DM. For light DM, m x < 30 GeV, the correct relic density 
requires a non-perturbatively large coupling g x so that the predictions should be taken 
only as 0(1) estimates in that region. Note that all these non-perturbative regions are, in 
addition, excluded by bounds on the decay width of the Higgs into invisible final states (see 
the discussion in Section 11C). 

In Figs. 3, 5, 7, and 9, we show the relic abundance curves for g x as a function of m x for 
the case where the light Yukawa couplings saturate their upper experimental bound. The left 
panels show the case where k u = 0.98 mb/m u and all the other couplings at their SM values, 
the middle panels the case where Kd = 0 .93rrib/md, and the right panels the case where 
k s = 0.70 m^/m s . In all of these cases the cross section for DM annihilation to light jets, 
a(xx jj)t coming from DM annihilating to light quarks, is comparable to the annihilation 




























































































Scalar DM, k u = O.Q&mb/m u 


m x [GeV] 


Scalar DM, k d = 0.93mj/mj 


m x [GeV] 


Scalar DM, n s = 0.70 mi,/m s 


m x [GeV] 


Figure 3: Bounds on the Higgs-portal coupling g x for scalar DM, assuming maximal allowed values 
for the Yukawa couplings to the u, d, s quarks (left to right), keeping all the other couplings to 
their SM values. The color coding is the same as in Fig. 2. 



Figure 4: Bounds on the Higgs-portal coupling for vector DM, assuming vanishing (left) and SM 
(right) Yukawa couplings to u, d, s quarks. The color coding is the same as in Fig. 2. 

cross section to 6-jets, cr(xX 66). For m x < mw these are the two dominant annihilation 
modes. Since the annihilation cross sections to 6-jets and light jets are comparable, the relic 
abundance curve show only a small change in g x when the 6-quark threshold is reached. This 
should be compared with the case of the SM Yukawa couplings shown in the right panels 
of Figs. 2, 4, 6, and 8. In this case the XX 33 annihilation is almost exclusively due to 
DM annihilating to two gluons, so that cr(xx 33 ) ^ a (xx 66), while the annihilation 
into two light quarks is negligible. For m x < mw and SM Yukawas, the required g x is thus 
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Vector DM, = 0.93mf,/m,j Vector DM, k s = 0.70m;,/m s 


1 10 10 2 10 3 1 10 10 2 10 3 10 4 
m x [GeV] m x [GeV] 


Vector DM, k u = 0.98 nn\,/m u 


10 

m x [GeV] 


Figure 5: Bounds on the Higgs-portal coupling for vector DM, assuming maximal allowed values 
for the Yukawa couplings to the u, d, s quarks (left to right), keeping all the other couplings to 
their SM values. The color coding is the same as in Fig. 2. 




m x [GeV] m x [GeV] 


Figure 6: Bounds on the Higgs-portal coupling for Dirac DM, assuming A = 1 TeV and vanishing 
(left) and SM (right) Yukawa couplings to u , d, s quarks. The color coding is the same as in Fig. 
2 . 

bigger by 30% — 40% than in the case of light Yukawas at their present experimental limits, 
and exhibits a significant jump below the 5-quark threshold. 
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Figure 7: Bounds on the Higgs-portal coupling for Dirac DM, assuming A = 1 TeV and maximal 
allowed values for the Yukawa couplings to the u, d, s quarks (left to right), keeping all the other 
couplings to their SM values. The color coding is the same as in Fig. 2. 




m x [GeV] m x [GeV] 


Figure 8: Bound on the pseudoscalar Higgs-portal coupling for Dirac DM, assuming A = 1 TeV 
and vanishing (left) and SM (right) Yukawa couplings to u, d, s quarks. The color coding is the 
same as in Fig. 2. 
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Dirac DM, CPV, k u = 0.98rrif,/m„ 


m x [GeV] 


Dirac DM, CPV, kj = 0.93mj,/mrf 


m x [GeV] 


Dirac DM, CPV, k s = 0.70 mh/m s 


m x [GeV] 


Figure 9: Bounds on the pseudoscalar Higgs-portal coupling for Dirac DM, assuming A = 1 TeV 
and maximal allowed values for the Yukawa couplings to the u, d, s quarks (left to right), keeping 
all the other couplings to their SM values. The color coding is the same as in Fig. 2. 

C. Invisible decay width of the Higgs 


The bounds on the invisible decay width of the Higgs boson provide stringent constraints 
on Higgs-portal DM [7]. The partial h —» XX decay widths are given by 


ri, = 


al '"tv 




xx 167r M h r x 

8iT 


r?. F =|M4^ + | 


XX 


(17) 


rL = 


9x M h v w 


xx 1287r m* 


Px (! - r x + ¥i) > 


where r x = 4 m^/Ml and /3 X = y/1 — r x . 

The current best limits on the invisible branching fraction of the SM Higgs are obtained 
from Zh production. The CMS collaboration gives a 95% CL limit of Br {h — y inv) < 0.58 for 
M h = 125 GeV [42] and ATLAS finds Br (h —> inv) < 0.75 for M h = 125.5 GeV [43]. Note 
that the increased light-quark Yukawa couplings, at their presently allowed values, do not 
appreciably change the Higgs production cross section [13]. Their main effect is to increase 
the total decay width of the Higgs and thus reduce the branching ratios to the other decay 
modes: 

r(ft -»■ xx) 


Br(h —> xx) = 


(18) 


r (h ->■ xx) + r^ ot x [l + - i)Br SM 0 -»■ qq)] ' 

In Figs. 2 to 9 we denote the bound on g x corresponding to the ATLAS upper limit on 
Br {h —y inv) with a dashed black line and grey out the excluded region in the g x vs. m x 
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Figure 10: The effect of large light-quark Yukawa couplings on the indirect detection bounds from 
Fermi-LAT observations of Milky Way dwarf spheroidal satellite galaxies [45] . 

plane. We see that the light DM Higgs portal, m x < m^/2 is excluded by the Higgs invisible 
decay width. 

Vector boson fusion, gluon fusion and ttH production, with the off-shell Higgs going to 
two DM particles, can provide some limited sensitivity to DM masses above rrih/2. A com¬ 
bination of the searches in the three channels at a 100 TeV collider could exclude the scalar 
thermal relic DM Higgs portal for DM masses in parts of the m^f 2 < m x < niw interval at 
95% C.L. [10] (these result receives only a negligible correction if light quark Yukawa cou¬ 
plings are enhanced). For m x < rrih /2 the invisible Higgs decay width is, however, always 
the most constraining [44], 

D. Indirect detection 

In indirect signals of DM annihilation, the effect of changing the light-quark Yukawa 
couplings within the presently experimentally allowed ranges leads to at most 0(1) effects. 
Further, the effect is present only for DM masses below the W threshold where the dominant 
annihilation channel is into the bb final state. For example, Fig. 10 shows the recast of 
the Fermi-LAT bound from dwarf spheroidals [45] for scalar DM, following the procedure 
outlined below. 

Photon flux measurements with y-ray telescopes can put a strong bound on the annihila¬ 
tion cross-section of the DM. The strongest bound in the DM mass range of interest has been 
recently released by the Fermi-LAT collaboration [45] based on Pass 8 observation data of 
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the Milky Way dwarf spheroidal satellite galaxies (dSphs). There is also a recent analysis 
based on the Dark Energy Survey (DES) dSph candidates using the Fermi-LAT data [46]. 
While this bound is competitive with the one from the known dSphs, it is still weaker on 
its own. One could also consider the bounds from the isotropic gamma ray background 
(IGRB) [47]. In our analysis, we recast the Fermi-LAT bound on the bb final state using a 
simple re-weighting procedure of the photon spectra which will be discussed below. 

The observed differential photon flux from the annihilation of dark matter is given by 


d<f> 


dErydQ 


1 1 

47t 2 m x 


J 


/ \ °^7 

2^( av )f 

L / 


dE 


(19) 


where J is an astrophysical factor which depends on the distance to the source and the 
dark matter density profile. The factor in the brackets is the one most interesting for our 
purposes. It depends on the velocity-averaged cross-section and photon spectrum per DM 
annihilation. 

The Fermi-LAT analysis gives bounds for the different final states separately while we 
have an admixture of final states. In order to recast the bound, we rely on the observation 
that for heavy DM the photon spectra from DM annihilation into quarks, gauge bosons, and 
the Higgs boson all peak at approximately the same photon energy and have approximately 
the same shape. Therefore, to extract the bound on the DM Higgs portal coupling g x , it is 
sufficient to find the zeros of the polynomial 


f(9x) = J2( av )f 

f 


dNf 

dE 


(av) 


FERMI 

b 



( 20 ) 


where in the last term, (ctu)™ rmi is the bound from the Fermi-LAT analysis on the velocity- 
averaged cross-section. The photon spectra were obtained from the interpolation tables 
provided in the PPPC4DMID package [48]. In all cases except for the hh final state, f(g x ) 
has only one zero up to a sign ambiguity. For the hh final state, however, the zero of f(g x ) 
closest to the g x corresponding to \X ~► bb is the one used to rescale the Fermi-LAT bound 
on g x as a function of m x . The resulting bounds are shown in Figs. 2 to 9. 


E. Direct detection 

We have seen so far that most DM observables exhibit only a weak dependence on the 
light-quark Yukawas. This is not the case for the direct DM detection. In fact, modifying 
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the light-quark Yukawa couplings can significantly change the predictions for DM - nucleus 
scattering cross sections. 

The differential cross section for spin-independent DM scattering on a nucleus is given 

by _ 

da = m A \M\ 2 
d,E R /J^t’rei 327TS ’ 

where E R is the nuclear recoil energy, tyia is the mass of the nucleus, h x a = m x mA/(m x -\-rriA) 
is the reduced mass of the DM - nucleus system, s = ( m x +mA ) 2 is the center-of-mass energy, 
u re i is the DM velocity in the detector rest frame, and \A4\ 2 is the spin-averaged squared 
matrix element. 

The matrix element M. depends on the effective Higgs couplings to the nucleus. Since 
the momentum exchanges in DM scattering on nuclei are much smaller than the Higgs mass, 
we can calculate \M. | 2 by first integrating out the Higgs and the heavy quarks (t, b , c). This 
gives an EFT with light quarks and gluons interacting with DM through local operators, 
described by the effective Lagrangians 



Eq — QQ-C, 




q - 


-G a ^G a ^ . 


Here, the last two terms arise from CP -violating Higgs couplings. The sums are over the 
light quarks q — u,d, s. The heavy quarks are integrated out and contribute only via the 
gluonic terms in the current. For the two corresponding dimensionless Wilson coefficients 
we have 

Cg = + Kt + Kb + H c , Cg = + Kt + Kb + R c , (26) 


where the first contribution is from tree-level matching, and the remaining from one-loop 
matching, working in the limit of heavy quarks. This is well justified for top and bottom 
quarks. For scattering on heavy nuclei, e.g., on Xe or W, the maximal momentum exchanges 
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for DM with mass above approximately 1 TeV may, however, start to become comparable 
to the charm-quark mass. We neglect these effects, while they may need to be included in 
the future if such heavy DM is discovered. 

CP-violating Higgs couplings to light quarks lead to spin-dependent interactions of DM 
with the target nuclei. The corresponding scattering rates are suppressed relative to the 
spin-independent interaction rates from CP-conserving Higgs couplings. We will therefore 
neglect the CP-violating interactions in our numerical analysis of direct detection scattering 
rates; i.e., we will set R q = 0, C q = 0 from now on. 

The nucleon matrix elements of the remaining terms in the scalar current S q are conven¬ 
tionally parametrized by (see, e.g., [49]), 

(N\m q qq\N) = m N f ^ N q ] , (27) 

(Nl-^G^anN) = —rUN^TC ■ (28) 

In the heavy-quark limit for t , b, c the trace anomaly equation leads to the relation [49, 50] 


= 1 - E C' (29) 

q=u,d,s 

We can also define the effective Higgs coupling to nucleon as the expectation value of the 
scalar current, 


/f > = (iV|S„|JV) = 


m N 
V\y L27 


Cg + ( K <? 27 ^ 9 ) f 


AN) 

'Tq 


(30) 


The exclusion curves from LUX [51] and XenonlOO [52], assuming a local DM density 
of 0.3 GeV/cm 3 , are shown in Figs. 2 to 9 as blue and red bands, respectively. The width 
of the exclusion curves represents the uncertainties in the hadronic matrix elements of the 
light-quark scalar currents. For the s quark we use the lattice determination = 0.043 ± 
0.011 [53]. The matrix elements for u and d quarks can be related to the ny-jv term. A 
Baryon Chiral Perturbation Theory (ByPT) analysis of the nN scattering data gives ny-jv = 
59(7) MeV [54], This is in agreement with ByPT fit to world lattice Nf = 2 + 1 QCD 
data, which gives oyjv = 52(3)(8) MeV [55]. Including both A(1232) and finite spacing 
parametrization in the fit shifts the central value to cr n N — 44MeV. To be conservative we 
use oyTv = (50 ± 15)MeV, which gives = (1.8 ± 0.5) x 10~ 2 , = (3.4 ± 1.1) x 10 -2 , 

= (1.6 ± 0.5) x 10~ 2 , f^l = (3.8 ±1.1) x 10 2 , using the expressions in [56, 57]. For the 
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effective Higgs coupling to nucleons this gives 


Ap) 

Js 


An) 

Js 


mw 

V\v 


m w 

v w 


(1.8 ± 0.5)k u + (3.4 ± l.l) Kd + (4.3 ± 1.1 )k* 


+ (6.70 ± 0.12) (k c + K b + n t + ^ p ) 


x 10 


-2 


(1.6 ± 0.5)k„ + (3.8 ± l.l)K d + (4.3 ± l.l)n s 


+ (6.69 ± 0.12) (k c + n b + n t + 


x 10 


-2 


( 31 ) 


(32) 


We use the results in [58] to relate the nuclear matrix elements to actual scattering rates on 
nuclei via nuclear form factors. 

We show the direct detection exclusion limits for SM (n u ^ s = 1) or vanishing (n U) d )S = 0) 
light-quark Yukawa couplings in the right and left panels in Figs. 2, 4, 6, and 8, respectively. 
The exclusion limits are approximately two times weaker in the latter case; the constraint 
does not vanish because, for small values of the light-quark Yukawas, the scattering cross 
section is dominated by the gluon part of the scalar current, Eq. (25). When the light- 
quark Yukawas are taken to be at the upper limit of their experimentally allowed range, 
i.e. comparable to the SM bottom Yukawa, the direct detection bounds on g x become 
significantly stronger, by a factor of about m b /m q (Figs. 3, 5, 7, and 9). 

It is interesting to note that, because of the dominance of the gluon contribution, for small 
light-quark Yukawas the theory uncertainty in the exclusion bands is significantly smaller 
than if the light Yukawa couplings are allowed to saturate the present experimental bounds. 
(The nuclear matrix element of the effective gluon term has smaller relative uncertainties 
than the corresponding matrix elements of m q qq since 

For m x smaller than a few TeV, the DM direct detection bounds are compatible with 
the thermal relic Higgs-portal DM only if light quark Yukawas are well below the present 
experimental bounds (the exception is a pseudoscalar fermion DM with enhanced strange 
Yukawa, where the bound is m x > m^/ 2, see Fig. 9). This means that if thermal relic DM is 
discovered, it would immediately place an upper bound on k u , hid, K a , assuming Higgs-portal 
mediation (unless in the case of fermion DM that has purely pseudoscalar couplings). We 
comment in more detail on that observation in Section IV. 

Since the DM - nucleus scattering cross section is the only DM observable that exhibits 
a rather pronounced dependence on the values of the light-quark Yukawas, we study this 
dependence in more detail. 
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Figure 11: Left: The ratio of direct detection bounds on g x from Xenon target varying k u (dark 
red), Kd (light red), or k s (blue), and the bound on g x assuming SM Higgs Yukawa couplings. The 
LHC upper bounds on /q are denoted by vertical dashed lines with shaded regions excluded. Right: 
the ratio of predicted scattering cross sections. The dotted lines correspond to negative values of 


In Fig. 11 (left) we show how the direct detection bounds on g x are affected by changes 
in the values of the light-quark Yukawas. We plot the ratio 


Gx — 


9y 


(33) 


5 9x 5f“ ax ( 1) ’ 

where (?™ ax (fi: 9 ) is the upper bound on g x obtained from direct detection experiments for 
a given value of K q , with q = u,d, s. Hence, g™ ax (l) is the bound obtained assuming SM 
Yukawa couplings. Its value depends on m x , on whether DM is a scalar, fermion, or vector, 
and on the experiment that measured the bounds. Similarly, also q x ax (K q ) depends on m x , 
the spin of DM, and the experiment; however, all these dependences cancel in the ratio £ 3x . 
The ratio £ 9x thus only depends on K q and on which target material was used to derive the 
direct detection bounds. In Fig. 11 we show £ 9x for a Xenon target, varying in turn k u (dark 
red line), Kd (light red) and k s (blue), while keeping all other parameters fixed to their SM 
values. We set the hadronic matrix element to their present central values, anticipating 
that in the future their uncertainties will be further reduced. In Fig. 11 (right) we show a 
closely related quantity - the ratio of the scattering cross sections with varied k u ^, s and the 
scattering cross section with SM Yukawa couplings, cr d . d ./cr|^ . 
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Fig. 11 illustrates clearly that the difference between the bounds where one of the light 
quark Yukawa couplings is taken to be small or vanish completely, and the bounds where 
all the couplings are SM-like, is very small, 0(10%). Saturating the present experimental 
bounds on k u or Kd, the allowed value of g x could lie two orders of magnitude below what 
one obtains for the case of SM Yukawa couplings. Such large values for the light-quark 
Yukawas are not very likely to be realized in a concrete model, as we will discuss in the 
next section. However, it is very interesting to observe that even a moderate increase of the 
values of the light-quark Yukawa couplings to only a few times their SM values can have a 
significant effect on the direct detection bounds, enhancing the scattering cross sections by 
up to a factor of ten. 

Finally, in the right panel of Fig. 11 we show the cross section ratios for negative values 
of the n q (dotted lines). We see that the interference of the light-quark contributions with 
the effective gluon interaction can, in principle, lead to a strong reduction of the scattering 
cross section. 

III. CHANGES TO YUKAWA COUPLINGS IN NEW PHYSICS MODELS 

So far we allowed the Yukawa couplings of the Higgs to quarks to have arbitrary values, 
only restricting them to lie within the bounds obtained from global fits to LHC data. For 
simplicity, we also neglected flavor violation and CP violation when discussing their impact 
on the DM interactions. 

Of course, changes of the Yukawa couplings by several orders of magnitude, as allowed 
by current experimental constraints on the light-quark Yukawas, are not very likely to be 
realized in a complete model, and might require significant fine tuning of the corresponding 
quark masses. In this section we investigate how large the deviations from the SM Yukawa 
interactions can be in popular models of NP with viable flavor structures. 

Tables I and II summarize the predictions for the effective Yukawa couplings in the Stan¬ 
dard Model (SM), in multi-Higgs-doublet models (MHDM) with natural flavor conservation 
(NFC) [59, 60], in the MSSM at tree level, the Giudice-Lebedev model of quark masses 
(GL) [61], in NP models with minimal flavor violation (MFV) [62], in Randall-Sundrum 
models (RS) [63], and in models with a composite Higgs, realized as a pseudo-Nambu- 
Goldstone boson (pNGB) [64-67]. For completeness, we include both the flavor-conserving 
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Table I: Predictions for the flavor diagonal up-type Yukawa couplings in a number of new physics 
models (see text for details). 


and flavor-violating Yukawa interactions, and allow for CP violation. The Higgs couplings 
to quarks are thus described by 


m Q- m q- i T / x _ . 

£cff,q = -Kq—qqh - iK q — q^qh - ( K qq ' + tn qq :)q L q R h + h.c. 
I’w Vw 1 


(34) 


where a sum over the SM quark fields is understood. The first two terms are flavor diagonal, 
with the first term CP conserving and the second term CP violating, and coincide with 
the definitions in eqs. (3) and (5), respectively. The terms in square brackets are flavor 
violating, with the real (imaginary) part of the coefficient CP conserving (violating). In the 
SM we have K q = 1, while k q = K qq i = R qq i = 0. The flavor-violating couplings in the above 
set of NP models are collected in Tables III and IV. These tables complement the analyses 
in [68-70] (see also [71], where implications of a negative top-quark Yukawa were explored, 
and [72] for an indirect bound on the down-quark Yukawa in alignment models). 


A. Dimension-Six Operators with Minimal Flavor Violation 

We start our discussion by considering dimension-six operators arising from integrating 
out NP at a high scale A. In addition, we assume that the flavor breaking in the NP sector 
is only due to the SM Yukawas, i.e. that NP satisfies the Minimal Flavor Violation (MFV) 
hypothesis [62, 73-78]. Integrating out the new physics states gives for the Higgs couplings 
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Table II: Predictions for the flavor diagonal down-type Yukawa couplings in a number of new 
physics models (see text for details). 
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Table III: Predictions for the flavor violating up-type Yukawa couplings in a number of new physics 
models (see text for details). In the SM, NFC and the tree-level MSSM the Higgs Yukawa couplings 
are flavor diagonal. The estimates of the CP-violating versions of the flavor-changing transitions, 
Kij/nt , are the same as the CP-conserving ones, apart from substituting “Im” for “Re” in the 
“MFV” row. 


to quarks 

£eft = Y u Q L H c u r + Y d Q L Hd R + ^Q L H c u R (H*H) + Q L Hd R (H ] H) + h.c., (35) 

where A is the scale of new physics and H c = icr 2 iL*. We identify the NP scales in the up- 
and down-quark sectors for simplicity. There are also modifications of quark kinetic terms 
through dimension-six derivative operators. These can be absorbed in (35) using equations 
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Table IV: Predictions for the flavor violating down-type Yukawa couplings in a number of new 
physics models (see text for details). In SM, NFC and tree level MSSM the Higgs Yukawa couplings 
are flavor diagonal. The estimates of the CP-violating versions of the flavor-changing transitions, 
/Kb, are the same as the CP-conserving ones, apart from substituting “Im” for “Re” in the 
“MFV” row. 


of motion [79]. The quark mass matrices and Yukawa couplings after EWSB are thus 


M = V ^(Y +Y’^-\ 
q y/2V q + r *2A*)' 


y q = Y q + 3 Y’ 


J w_ 
q 2A 2 


q = u,d. 


(36) 


Because Y q and Y q > appear in two different combinations in M q and y q , the two, in general, 
cannot be made diagonal in the same basis. 

In MFV the coefficients of the dimension-six operators can be expanded in terms of Y u ^ 


(37) 


K = o»n + b u r u YjY u + c u Y d yJy u + ■■■, 
y; = a d Y d + b d Y d Yj Y d + c d Y u Y<Y d +■■■. 

with a q ,b q ,c q ~ C(l). Working to hrst order in dimension-six operator insertions we can thus 
write for the Yukawa couplings, in the mass eigenbases for up and down quarks respectively, 


Du 

Vd = 


1 + y 1 u + buiusu) 2 + c uV(y d SM ) 2 ^ H - 

i + (°cf + bd(y d SM ) 2 + c dV\ysM) 2 V + • • • 


VsM 


VSM 


(38) 


Here y^ are the diagonal matrices of the SM Yukawa couplings, while V is the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. In general, the coefficients a q , b q , c q are complex so that 
CP-violating Higgs couplings arise at 0{vyy/ A 2 ). Flavor-violating Higgs couplings arise first 
from the terms proportional to c U)< 2 in the Yukawa expansion and are thus suppressed by 
the corresponding CKM matrix elements. In Tables I-IV we collect the values for flavor- 
conserving and flavor-violating Yukawa couplings in the “MFV” row, assuming that all the 
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coefficients a q , b q , c q are 0( 1), and show only the numerically leading non-SM contributions. 
In the expressions we also set V tb to unity. 

The corrections to DM phenomenology are dominated by changes of the third-generation 
Yukawa couplings. The MFV corrections to light-quark Yukawa couplings are all either 
additionally CKM suppressed or involve extra insertions of light-quark masses. Hence the 
theory error in Higgs-portal DM phenomenology due to Yukawa coupling uncertainties will 
be small in MFV models of NP once the Higgs couplings to top and bottom quarks are well 
measured. 

B. Multi-Higgs-doublet model with natural flavor conservation 

In MHDMs there are no tree-level FCNCs if natural flavor conservation is assumed [59, 
60]. Under this assumption we can choose a Higgs doublet basis in which only one doublet, 
H u , couples to the up-type quarks, and only one Higgs doublet, H d , couples to the down- 
type quarks 2 . After EWSB the two doublets obtain the vevs v u and v dl respectively. On the 
other hand, the vevs of all Higgs doublets contribute to the W and Z masses. They satisfy 
the sum rule v w = Ei v h where the sum is over all Higgs doublets. 

The neutral scalar components of Hi are (vi + hi)/y/2, where the dynamical fields hi are a 
linear combination of the neutral Higgs mass eigenstates (and include h u and h d ). We thus 
have hi = V bl h + ..., where V bl are elements of the unitary matrix V that diagonalizes the 
neutral-Higgs mass terms and we only write down the contribution of the lightest Higgs, h. 
Under the assumptions above, the mass and Yukawa terms can be diagonalized in the same 
basis, so that there is no flavor violation and no CP violation in the Yukawa interactions: 

K qq ! R qq ' 0, R q 0. (39) 

We obtain a universal shift in all up-quark Yukawa couplings, and a different universal shift 
in all down-quark Yukawa couplings, given by 

K u = K c = K t = V hu — , K d = K s = K b = V hd — . (40) 

Vu V d 

Since the shifts are universal over generations, the precise measurements of the Higgs cou¬ 
plings to top and bottom quarks will also determine the Higgs couplings to light quarks. 

2 Note that H u = H d is included as a special case. 
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Both Kt and Kb are expected to be known with 0(5%) precision after the end of the high- 
luminosity LHC run [80, 81]. The uncertainties in the DM direct detection rates due to 
uncertainties in the Yukawa couplings will thus be negligible, assuming NFC. Note that the 
Higgs portal with an additional SM singlet mixing with the Higgs is also described by the 
above modifications of fermion couplings, with a completely universal shift Ki = cos 0, where 
0 is the singlet-Higgs mixing angle [82, 83]. 

Our analysis of modified Higgs-portal DM phenomenology given in Section II applies in 
the somewhat special limit where the DM only couples to the lightest mass-eigenstate h. 
For instance, for scalar DM the general Higgs portal is 

£nfc = 9x,ijX^xH}Hj. (41) 

If the hcrmitian matrix of couplings g Xjl] is such that it has h as the only eigenstate with 
nonzero eigenvalue, then our analysis in Section II applies unchanged. In general, however, 
all the expressions in Section II get corrected by terms of order l/m 2 H . due to exchanges of 
heavy Higgs bosons with masses m^. If DM is heavy, m x > rriHi, the presence of heavy 
Higgs bosons would also open new annihilation channels. 


C. Type-II Two-Higgs-Doublet Model 


The MSSM tree-level Higgs potential and the couplings to quarks are the same as in the 
type-II two-Higgs-doublet model (2HDM), see, e.g., [84], This is an example of a 2HDM 
with natural flavor conservation in which v u = sin /3vw, Vd = cos /3vw- The mixing of h u j, 
into the Higgs mass-eigenstates h and H is given by 



where h is the observed SM-like Higgs. Thus 


Hid 


COS O' 

sin/T 
sin a 
cos/3’ 


(42) 


(43) 


while the flavor-violating and/or CP-violating Yukawas are zero. In the decoupling limit 
(/3 — a = 7t/ 2) the heavy Higgs bosons become infinitely heavy, while the Yukawa couplings 


24 




tend toward their SM value, K t — 1. The global fits to Higgs data in type-II 2HDM already 
constrain /3 — a ~ 7t/ 2 [37, 85, 86] so that in this case the corrections to Higgs-portal DM 
phenomenology due to non-standard Higgs Yukawa couplings are small. 

As in the case of MHDM, the DM phenomenology of Section 11 remains unchanged only 
in the limit where the DM couples to the light Higgs h but not to the heavy Higgs H. In the 
general case, our analysis gets corrections that are relatively suppressed by 0(m\/m 2 H ). If 
we are not too far away from the decoupling limit these corrections can be neglected, while 
in parts of the parameter space, where cancellation can occur, the extra contributions are 
numerically important [8]. 


D. Higgs-dependent Yukawa Couplings 


In the model of quark masses introduced by Giudice and Ledebev (GL) [61] the Higgs- 
quark interactions are written in terms of effective operators 


£„ = c) 


(H*H\ 


fWH\ 


\ M 2 ) Q L > iUR jH + c ij ^ jyj - 2 J Ql, idRjH + h.c 


(44) 


They can be thought of as arising from integrating out heavy mediators at a large mass 
scale M. In this model the light quarks couple to the Higgs only through operators with 
mass dimension higher than four, i.e., for light quarks we have n- - d ^ 0. The values of the 
integers n“’ d , and of the coefficients that take values of order unity, are chosen such that 
the hierarchies of the observed quark masses and mixing angles are explained, after EWSB, 
in terms of the expansion parameter e = v^y/M 2 & 1/60. Thus, the Yukawa couplings are 
of the form 

I/SM + l)(ylf)sM. (45) 


After mass diagonalization the SM Yukawas are diagonal in the same basis as the quark 
masses, {y“- d ) sm oc 6ijm%’ d , while the y^ d are not diagonal in the same basis 3 . Using the 
ansatz = a; + bj’ d with a = (1,1,0), b d = (2,1,1), and b u = (2,0,0), this gives the 
deviations in the Yukawa couplings collected in Tables I-IV in the row denoted by “GL”. 
Since the couplings to the bottom quark is enhanced by a large factor, n b ~ 3, the simplest 


3 Note that the mixing of contributions from different effective operators that may have large relative phases 
could lead to sizeable CP-violating contributions to the Yukawa couplings. 
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version of the GL model is already excluded by the Higgs data on h —> WW, h —> ZZ and 
h — > 77 decays. 

We therefore modify the initial GL proposal and assume that we have two Higgs doublets 
in (44), H u that only gives masses to up-type quarks and Hd that only gives masses to 
down-type quarks. The correct mass and CKM angle hierarchy is obtained by using b d = 
(1, 0, 0) in the ansatz for n d g , and leaving a and b u unchanged. This gives satisfactory Higgs 
phenomenology at present as long as Kf, = sin a/ cos /3 ~ 1 up to (9(20%). In this limit 
also K t = cos a/ sin/3 ~ 1. The scaling of Yukawa couplings for this modification of the GL 
model is shown in Tables I-IV in the row denoted by “GL2”. 

In the GL model it is natural that the Higgs is the only state that couples to DM. The 
GL model is thus an example of Higgs-portal DM where the light-quark Yukawa couplings 
can substantially differ from their SM values. For instance, in GL2 k u ~ 7/sy, Hd — 5 Kb, 
k s ~ 3 Kb, k c ~ 3 K t . The coupling of DM to gluons (25) C g ~ 4 K t + Kb, so that C g ~ (5/3)C^ M , 
and C g ~ 0{C g ). Taking Kb ~ 1, this means that the effective Higgs coupling to nucleons, 
governing the direct DM detections rates, gets enhanced compared to the SM Higgs Yukawa 
couplings by 

f(p) t( n ) 

~4t —~ 1.2*7 + 1.3K b ~ 2.5, -44 —~ 1.3«t + 1.3/Cft ~ 2.6. (46) 

f(p) | A n > | v ' 

Js Ism Js Ism 

Here most of the enhancement over the SM comes from enhanced k u and Kd, which is also 
the reason for enlarged isospin breaking (the difference between and f^)- As a result of 
larger couplings to light quarks the spin-independent DM scattering cross section can thus 
be enhanced by an order of magnitude in the GL2 model of light-quark masses. 

E. Randall-Sundrum models 

The Randall-Sundrum (RS) warped extra-dimensional models with the SM fields propa¬ 
gating in the bulk provide a solution to the hierarchy problem and simultaneously explain the 
hierarchy of the SM fermion masses without large hierarchies in the initial five-dimensional 
(5D) Lagrangian [63, 87, 88]. The fermion zero modes are either localized toward the UV 
brane (for lighter fermions) or toward the IR brane (the top, the left-handed b quark and 
potentially the right-handed c quark) [89, 90]. The Higgs held and the Higgs vev are both 
localized toward the IR brane. Integrating out the Kaluza-Klein (KK) modes and working 
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in the limit of a brane-localized Higgs, the SM quark mass matrices are given, to leading 
order in v w! m 2 KK , by [91] (see also [92-100], and Ref. [101] for a bulk Higgs scenario) 


M d[u) 

*? 


[ F q Y l( 2 ) F d(u)\ i:j Vw ■ 


(47) 


Here, ttikk is the KK mass scale. The F q ^ Ut d are diagonal 3x3 matrices of fermion wave- 
functions for the left-handed electroweak quark doublets and the right-handed electroweak 
up and down quark singlets, respectively, evaluated at the IR brane. Assuming flavor an¬ 
archy, the 5D Yukawa matrices for up and down quarks, Yf 5 ^, are general 3x3 complex 
matrices with 0(1) entries. For a Higgs held propagating in the bulk, 5D gauge invariance 
guarantees Y[ iD = Y 2 5jD [91]. 

At leading order in v w/ m 2 KK the Higgs Yukawas are aligned with the quark masses, i.e., 


M uA = y u ,d + 0(v 2 w lm 2 KK ). (48) 

The misalignment arises from dimension-six operators that are generated by tree-level KK 
quark exchanges. They give 



where Y is a typical value of the dimensionless 5D Yukawa coupling and is in numerical 
analyses typically taken to be below Y < 4 (see, e.g., [95]). The Higgs mediated FCNCs are 
thus suppressed by the same zero-mode wave-function overlaps that also suppress the quark 
masses, giving rise to the RS GIM mechanism [102-104], 

Using that the CKM matrix elements are given by V l0 ~ F qi /F qj for i < j, Eq. (49) can 
be rewritten as 


[Uu{d)\ [M u>d \ t 


--Y- 


J w 


l KK 


"Sigj) V 

v w V V ’ 


YiiMw- 1 

v w ij ’ 


* < J , 
i = 3, 
j <i. 


(50) 


This yields the Ki collected in Tables I-IV. In the numerical analysis of ref. [91] the diagonal 
values Kj were typically found to be smaller than one, with deviations in K t up to 30%, Kb 
up to 15%, in k S)C up to ~ 5%, and in k u ^ of 1% (these estimates were obtained fixing 
the mass of the Erst KK gluon excitation to 3.7 TeV, above the present ATLAS bound 
[105]). The effective Higgs coupling to nucleons, f$ N \ thus only gets reduced by 0(10%), 
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giving a 0(20%) smaller DM scattering cross section on nuclei, compared to the case of SM 
Yukawa couplings. The largest effect arises in DM annihilations to top quarks, where the 
cross section can be reduced by a factor of two, while the annihilation cross section to bb 
pairs can be ~ 30% smaller than for SM Yukawa couplings. 

F. Composite pseudo-Goldstone Higgs 

Finally, we investigate the possibility that the Higgs is a pseudo-Goldstone boson arising 
from the spontaneous breaking of a global symmetry in a strongly coupled sector [64-67]. 
We assume that the SM fermions couple linearly to composite operators Ol,r [106], 



(51) 


where i,j are flavor indices. This is the 4D dual of the fermion mass generation in 5D RS 
models. The Higgs couples to the composite sector with a typical coupling y*. The SM 
masses and Yukawa couplings then arise from expanding the two-point functions of the Ol,r 
operators in powers of the Higgs held [107], giving rise to four- and higher-dimensional Higgs 
operators, such as in (35). 

The new ingredient, related to the pNGB nature of the Higgs, is that the shift sym¬ 
metry dictates the form of the higher-dimensional operators. The flavor structure and the 
composite Higgs coset structure completely factorize if the SM fields couple to only one 
composite operator. The general decomposition of Higgs couplings then becomes [107] (see 
also [108, 109]) 


y,/),.//«„ + v ,; q , v , - ... -> pih/f )«//«;, , 


(52) 


and similarly for the down quarks. Here P(h/f) = ao + « 2 (h//) 2 + ... is an analytic function 
whose form is fixed by the structure of the spontaneous breaking and the embedding of the 
SM fields in the global symmetry of the strongly coupled sector, while / is the equivalent 
of the pion decay constant and is of order v\y Since the flavor structure of the coefficients 
of the dimension-four and dimension-six operators is the same, they can be diagonalized 
simultaneously. All corrections to the quark Yukawa couplings from this effect are therefore 
strictly diagonal, and we have 
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( 53 ) 



For example, for the models based on the breaking of SO( 5) to 5'0(4), the diagonal Yukawa 
couplings can be written as [110] 

1 + 2m - (1 + 2m + n)(v w / f) 2 


K q = 


(54) 


\J 1 ~ i v w/f) 2 

where n,m are positive integers. The MCHM4 model corresponds to m = n = 0, while 
MCHM5 is given by m = 0, n = 1. 

The flavor-violating contributions to the quark Yukawa couplings then arise only from 
corrections to the quark kinetic terms. That is, they are related to dimension-six operators 
of the form [107] 

/yt /y /yt /y 

(55) 


- H*H _ H'H 
(IlW^L 7 9 , UrW U R~ 


A 2 ’ " ^ « A 2 , • ■ 

These operators arise from the exchange of composite vector resonances with typical mass 
M* ~ A. After using the equations of motion they contribute to the misalignment between 
the fermion masses and the corresponding Yukawa couplings. The NDA estimates for these 
corrections are, neglecting relative 0( 1) contributions in the sum [15, 107, 111], 


K,. . nu 
IJ 


2 


mu i , X u ,« m * 
- r An .A, 


V W 


'R,i / 'R,j 


V W 


(56) 


and similarly for the down quarks. If the strong sector is CP violating, then ~ 

The exchange of composite vector resonances contributes also to the flavor diagonal 
Yukawa couplings, shifting the estimate (53) by (note the different normalizations of R q 
and K qq i in (34)) 


2 r 


Ar q . ~ 

Hi 


2y 2 — 
y *M 2 


W,T + 0«,<) : 


(57) 


This shift can be large for the quarks with a large composite component if the Higgs is 
strongly coupled to the vector resonances, y* ~ 47T, and these resonances are relatively 
light, M* ~ Attvw ~ 3 TeV. The left-handed top and bottom, as well as the right-handed 
top, are expected to be composite, explaining the large top mass (i.e., A| 3 ~ ^r ,3 ~ !)• 
In the anarchic flavor scenario, one expects the remaining quarks to be mostly elementary 
(so the remaining A* 1). However, if there is some underlying flavor alignment, it is also 
possible that the light quarks are composite. This is most easily achieved in the right-handed 
sector [109, 112, 113]. Taking all right-handed up-type quarks fully composite, and assuming 
that this results in a shift Ak u ~ Ar c ~ Ati t ~ 1, this would lead to an increase in the 
effective Higgs coupling to nucleons, fg, of about 50%, and an increase in the DM-nucleon 
scattering rate of about 100%. 
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IV. CONSTRAINING THE LIGHT-QUARK YUKAWA COUPLINGS 


If DM is a thermal relic interacting with ordinary matter predominantly via SM Higgs 
exchange, direct detection scattering rates immediately give information about the light- 
quark Yukawa couplings once the coupling of the DM particle to the Higgs particle is fixed. 

In fact, DM scattering in direct detection searches would be one of the very few possible 
probes of the light-quark Yukawa couplings. The interactions of the Higgs boson with u, d, 
or s quarks give rise to flavor-conserving neutral currents. Off-shell Higgs contributions in 
processes with only SM external particles always compete with other, much larger flavor- 
conserving neutral currents induced by gluon, photon, or Z exchange. This leaves us with 
two options: either to consider on-shell Higgs decays [13-15], or to use new probes, such as 
DM scattering in direct detection experiments. 

In principle, there is enough information to make a closed argument. Suppose that in¬ 
direct DM searches yield a positive DM annihilation signal for m x > rrih/Z. At the end of 
the high-luminosity LHC run, the Higgs couplings to W, Z, t, and b will be precisely deter¬ 
mined. Assuming that DM is a thermal relic interacting only through the Higgs portal, this 
fixes the value of g x since the annihilation cross section for m x > nth /2 is otherwise almost 
completely controlled by the Higgs couplings to W, Z, and t. In principle, a consistency 
check that the DM is really interacting through a Higgs portal could be provided, for a very 
limited range of DM masses m x > m^/2, by a 100-TeV hadron collider [10]. 

After the discovery of DM, the direct detection searches would immediately imply an 
upper bound on the light-quark Yukawa couplings. As an illustration, consider the excess 
in y-ray emmission in the recently discovered dwarf spheroidal galaxy Reticulum 2 [46]. 
Let’s take the bold step of interpreting this signal as originating from DM annihilating into 
bb pairs (see Ref. [114] for details). Assuming the Dirac-fcrmion DM scenario with purely 
CP-violating couplings, we obtain a la region in the m x - g x plane that is not yet excluded 
by direct detection constraints, denoted by the orange lines in Fig. 12. Note that part 
of this region is consistent with DM furnishing the dominant component of the observed 
relic density while at the same time not being excluded by the invisible Higgs decay width. 
Concentrating on the overlap region, m x ~ 75 GeV, a comparison with the ratios shown in 
Fig. 11 would immediately imply an upper bound of k u < 10, Kd < 10, n s < 12 from the 
LUX direct detection search (allowing only one of the Yukawa couplings to float at a time). 
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Figure 12: The 7 -ray excess in the recently discovered dwarf spheroidal galaxy Reticulum 2, inter¬ 
preted as a signal of DM annihilating into bb pairs, is shown as the black la contour (see Ref. [114] 
for details). The orange lines show the 95% CL exclusion limits at the 14-TeV LHC (solid line) 
and a prospective 100-TeV hadron collider (dashed line), obtained by rescaling the bounds given 
in Ref. [10]. The remaining color coding is the same as in Fig. 2. See text for more details. 

These estimates could potentially be loosened by uncertainties in the DM velocity profile 
and the local DM density. On the other hand, if DM is discovered in direct detection the 
relative size of the light-quark Yukawas could be probed by comparing scattering rates on 
different target materials. 

An additional cross check of our scenario could be provided by searches for DM production 
at hadron colliders. In Fig. 12 we denote the 95% CL exclusion limits, assuming 3000/fb 
of data, at the 14-TeV LHC by a solid orange line and at a prospective 100-TeV hadron 
collider by a dashed orange line. These curves have been obtained by converting the bounds 
in Ref. [10] to the case of Dirac DM using FeynRules [115] and MadGraph5 [116]. We see 
that, while the LHC will be sensitive to a part of the interesting region in parameter space, 
the scenario of m x = 75 GeV DM can be excluded only at a 100-TeV collider. 

V. CONCLUSIONS 

Not much is known experimentally about the couplings of the Higgs to light quarks. It is 
entirely possible that the Higgs couples only to the third generation of fermions. Experimen- 
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tally equally viable is the possibility that the light-quark Yukawas are significantly enhanced, 
up to (9(50) for k s and up to 0( 10 3 ) for k u and k ( j. Such extremely large enhancements are 
not natural from a model-building point of view as they require a large fine tuning of the 
light-quark masses, but at present cannot be excluded experimentally. 

Modified Yukawa couplings to light quarks could have implications for DM searches. In 
this paper we focused on Higgs-portal DM. We considered constraints on scalar, vector, and 
fermionic Higgs-portal models of DM from relic density, direct and indirect detection, and 
the invisible Higgs width. A central result of our analysis is that, for phenomenologically 
viable Higgs-portal DM, there is a relatively small change in the predictions between the 
case where the Higgs is SM like and the case where the Higgs couples only to the third 
generation of fermions. For direct detection this is a consequence of the fact that, for 
very small light-quark Yukawas, the scattering cross section is dominated by the effective 
Higgs-gluon coupling, which is obtained by integrating out the heavy quarks. For instance, 
setting all the light quark Yukawa couplings to zero reduces the direct detection scattering 
cross section by ~ 50% compared to the case where the light quark Yukawa couplings are 
assumed to have SM values. Similarly, the relic abundance and indirect detection signals 
are dominated by the heaviest kinematically open annihilation channels, diminishing the 
importance of Higgs couplings to light quarks. 

On the other hand, saturating the experimentally allowed values for the light-quark 
Yukawas, the DM direct detection rates can increase by four orders of magnitude compared 
to the case where the light-quark Yukawa couplings are kept at their SM values. Moreover, 
negative values of the light-quark Yukawa couplings can result in a strong reduction of the 
nucleon scattering rates. The changes in DM annihilation rates are much smaller. The anni¬ 
hilation of DM into light quarks is a subleading effect, unless m x < mw■ Even in this case, 
the dominant annihilation channel is into bb pairs, while the annihilation to light quarks 
can constitute at most an (9(1) fraction if the current experimental upper bounds on the 
light-quark Yukawa couplings are saturated. A Higgs-portal for DM in this mass range is 
excluded either by bounds on the invisible Higgs decay width or by indirect DM searches. 

We also investigated the expected sizes of corrections to DM phenomenology due to 
changes in Yukawa couplings in a number of new physics models. The largest deviation in 
expected DM scattering rate on nucleons was found for a modified Giudice-Lebedev model 
of light-quark masses where up to an order of magnitude enhancement due to corrections 
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to light-quark Yukawa couplings are possible. Similarly, an 0( 1) change of the scattering 
rate is anticipated in a pseudo-Goldstone Higgs scenario with composite right-handed light 
quarks while in RS models with anarchic flavor a reduction of about 20% can be expected. 
The effects in MFV models, multi-Higgs models with natural flavor conservation, and the 
type-II two-Higgs-doublet model (i.e., the tree-level Higgs sector of the MSSM), on the other 
hand, are expected to be much smaller. 

Finally, we point out that a discovery of Higgs-portal DM in indirect searches would 
immediately imply an upper bound on the light-quark Yukawa couplings due to the upper 
bounds in direct DM searches. 
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